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I ^1 ' Abstract 

Kinetics of model catalytic processes proceeding on inhomogeneous 
^ ■ surfaces is studied. We employ an extended mean-field model that 

. takes into account surface inhomogeneities. The influence of surface 

diffusion of adsorbent on the kinetics of the catalytic process is in- 
Tij" ' vestigated. It is shown that diffusion is responsible for differences in 

■ the reaction rate of systems with different arrangements of active sites. 
The presence of cooperative effects between inactive and active sites is 
demonstrated and the conditions when these effects are important are 
discussed. We show that basic catalytic phenomena on nonuniform 
surfaces can be studied with mean-field modeling methods. 
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1 Introduction 

Catalytic processes on spatially inhomogeneous surfaces are examples of het- 
erogeneous reacting systems where the heterogeneity means not only that 



^Presently at Fritz Haber Institute for Molecular Dynamics, Hebrew University, 
Jerusalem 91904, Israel; email: cwiklik@gmail.com 



1 



a reaction takes place on the interface but also that the catalyst surface 
is heterogeneous. Current advances in nanotechnology made it possible to 
produce nanopatterned surfaces with assumed geometric properties and to 
study their reactivity. In the case of these surfaces the heterogeneity oc- 
curs in the nanometer scale and in the presence of finite surface diffusion 
may strongly influence the kinetics of surface processes and hence play a key 
role in the surface reactivity [HE]- A proper description of the kinetics on 
nanostructured surfaces is, in most cases, impossible using only the standard 
phenomenological chemical rate equations [4j, l5j. The classical description 
assumes an ideal mixing of adsorbate particles and excludes the presence 
of both the concentration gradients and the spatial nonuniformities of the 
surface. To overcome these limitations two groups of methods are usually 
employed in theoretical studies: kinetic Monte Carlo simulations (KMC) 
where adsorbate particles can directly undergo all studied surface processes 
and extended mean-field modeling (MF) where spatial inhomogeneities are 
taken into account indirectly. In most of the former studies surface diffu- 
sion of adsorbate was either neglected or considered only as a rapid process. 
Even for uniform catalysts surfaces diffusion was neglected in early works, 
with the well-known Ziff-Gulari-Barshad model as an example [6]. In the 
case on nonuniform surfaces in early works, representing both MC and MF 
treatment, surface diffusion was not included exphcitly [2] or even completely 
neglected [3 [HI E] • The kinetics was also studied in the limit of rapid diffusion 
where local nonuniformities were not included [TUl [HI 1121 [131 [IS]- In recent 
works concerning nonuniform surfaces it is usual to include surface diffusion 
in MC treatment |3l [T71 [181 EH] whereas in MF models only approximate 
descriptions were employed (diffusion between two regions - see |20]) or a 
process with coverage dependent diffusion coefficient). However, these MF 
models were dealing only with uniform surfaces [T7|. Evans and co-workers 
investigated the role of finite surface diffusion on the poisoning of catalyst 
surfaces using MF methods but also for uniform surfaces [lll[T5]. In this com- 
munication an extended mean-field description based on the Master equation 
formalism and describing the kinetics of catalytic processes (including finite 
surface diffusion) on heterogeneous surface is presented. This formalism takes 
into account both a local arrangement of reactive and non-reactive surface 
sites and local arrangement of adsorbate. Therefore, it lets to include and 
investigate an influence of geometric surface heterogeneities on the kinetics 
of considered catalytic processes in the regime of finite surface diffusion. 
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2 Mean-field model 



The mean-field model introduced here was derived by averaging Master equa- 
tions describing a time evolution of the density of the probability of occur- 
ring different local adsorbate/surface configurations on the catalyst surface. 
A similar method was proposed previously by other authors, however, they 
did not study in detail the role of finite surface diffusion for non-uniform 
surfaces [2], [71 [21]. We employed a coarse grain model of a surface with a sur- 
face represented as a two dimensional lattice. In order to study non-uniform 
surfaces we considered 2 types of lattice sites: active and inactive in the 
surface reaction. Adsorption, desorption and diffusion were allowed to pro- 
ceed at each site, and reaction only on the reactive ones. We were studying 
surfaces with the increasing ordering in the active sites distribution, starting 
from the random distribution, then active sites were arranged in the form of 
stripes of increasing width (see Fig. [1]). Two model catalytic reactions were 
considered: unimolecular and bimolecular (the latter according to Langmuir- 
Hinshelwood mechanism). First we investigated the catalytic process with 
an unimolecular reaction of the type A ^ P. At a surface with two types of 
centers the process can be presented as: 



where: A(^g-^ - reactant particle in the gas phase, o - surface center inactive in 
the surface reaction, A^ads) ^ reactant particle adsorbed on an inactive surface 
center, * - active surface center. A* - reactant particle adsorbed on an ac- 
tive center, P(g) - product particle desorbed to the gas phase. Eqs. ([1]) and 
([2]) describe reversible adsorption of reactant on, correspondingly, inactive 
and active surface centers. Eq. ([3]) corresponds to surface diffusion of reac- 
tant particles between active and inactive surface centers. Eq. (jlj) depicts 
the unimolecular reaction proceeding on active centers and the immediate 
desorption of product to the gas phase (the assumption of immediate des- 
orption of product particles is widely used in modeling of surface reactions, 
for example, in ZifT-Gulari-Barshad model; generally, product particles may 
stay or diffuse on the surface and poison the catalyst, however, for many 
practically important catalysts, product desorption is much more efficient 



A^g) + * 



A* 




A* 

A* + o 




(1) 
(2) 
(3) 
(4) 
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than other surface processes and can be considered as a rapid phenomenon). 
This model includes information about the surface heterogeneity but only in 
the form of spatially averaged concentrations of reactant on both active and 
inactive centers. 

In order to describe the local structure of the surface we considered prob- 
abilities of occurring different one-site configurations: [o] - probability of 
occuring an inactive unoccupied site, [*] - active unoccupied, [Ao] - inactive 
occupied and [A*] - active occupied. We took into account all processes of 
creation and annihilation of these configuration including surface diffusion 
(for instance, an inactive unoccupied site can be generated either by des- 
orption of reactant particle from an occupied inactive site or by diffusion of 
a particle from an occupied inactive site onto an unoccupied site neighbor- 
ing to the considered one). The following set of equations describes a time 
evolution of one-site probabilities: 

Pdes,inact[A°] + Pdiff,IA[A O | *] - 

-Pads,inact[°] - Pdiff,Al[^ * l°] i^) 
Pdes,act[^*] + 

-Pads,act [*] - PdiffJA [A O \ *] (6) 

Pads,act[*] + Pdiff,IA[A O \*] - 

-Pdes,act[A*] - PdiffMA * |o] - Pr[A*] (7) 
Pads,inact[°] + Pdiff,Al[A * |o] - 

-Pdes,vnact[Ao\ - PdiffJA[A O |*] (8) 

where: Pij - probabilities (per time unit) of desorption/adsorption of reac- 
tant particle on the inactive/ active sites, Pdiff,Ai {Pdiff,iA) - probability of 
diffusion of A particle from an active (inactive) onto an inactive (active) site, 
Pr - probability of surface reaction. On right-hand sides of above equations 
there are probabilities of occurring of two-center lattice configurations which 
are equal to the probabilities of finding in the system the pairs of neighboring 
sites in the given states. For example, the term [Ao |*] depicts the probability 
of occurring the local configuration where an inactive occupied lattice site 
has an active unoccupied site in the nearest neighborhood. The changes of 
two-center terms which are present on the right-hand side of these equations 



dt 

d[*] 
~dt 

d[A*\ 
dt 

d[Ao\ 
dt 
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can be written as a function of three-center probabilities: 



dm 



f{[k\l],[m\n\o]) 



(9) 



dt 



The simplest way to solve such a hierarchy of equations is to use one-center 
approximation, where ~ [ijij]- However, this one-center description 
omits the presence of spatial correlations. It can be noticed that two-center 
configurations characterized by probabilities [y4*|o] and [y4o|*] may be 
found on the surface only if active and inactive sites neighbor each other 
and the presence of such configurations depends on the spatial arrangement 
of active centers. For example, in the case of randomly placed active cen- 
ters with a low value of 6act (ratio of active to inactive sites) almost each 
active site on the lattice has inactive neighbors. On the other hand, if spa- 
tial correlations in the sites arrangement are present, a significant fraction of 
active sites may have no inactive neighbors. Therefore, spatial correlations 
in the arrangement of active centers can be included by introduction of an 
'inhomogeneity parameter' to approximate two-center probabilities: 



where the inhomogeneity parameter w is an average probability that an ac- 
tive center neighbors an inactive one in the considered system (it depends 
on the number of active-inactive neighbors in the system) and w can be cal- 
culated for each arrangement of active centers. In this approximation only 
one-center concentration terms are explicitly present, however, information 
about two-center correlations in the arrangement of surface sites is also taken 
into account with the inhomogeneity parameter. Similar approximation was 
proposed by Jansen and Hermse [2]. We obtained an analytical solution for 
stationary states of eqs. ([5]) - ([8]) with the approximation given by eqs. (fTO!) 
and (fTTj) taking an additional assumption that the probability of diffusion 
between active and inactive regions is equal, i.e., Pdiff,Ai = PdiffjA (for 
analytical solutions see Appendices). 

We also studied the kinetics of bimolecular surface reaction of the type: 
2A —>■ P proceeding according to Langmuir-Hinshelwood scheme. In this 
case eq. has the form: 2 A* Pf^g) + 2* and two-center terms [A * \A*\ 
appear. We approximated them introducing an inhomogeneity parameter Wr 
which corresponds to the probability of occurring at least one active center 



[A * |o] 
[A o I*] 



(10) 
(11) 
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in the nearest vicinity of considered active site. Assuming, as previously, 
Pdiff,Ai = PdiffjA we obtained analytical solution for stationary states (see 
Appendices). 

3 Results 

During this study we focused on the influence of both the arrangement of 
surface sites and surface diffusion on the kinetics of catalytic process. There- 
fore, we assumed Pads,act = Pads,inact = Pads and Pdes,act = Pdes,inact = Pdes 

(i.e., all surface sites had the same adsorption/desorption properties) and 
Pdiff,Ai = PdiffjA = Pdiff (due to the conditions for our analytical solu- 
tions). In order to test the mean-field approach presented above we com- 
pared mean-field results with the results obtained in our earlier work for 
similar system where Monte Carlo simulations were used [22]. Both surface 
coverages and turn-over numbers (results not presented here) were in a very 
good quantitative agreement with MC results in the case of unimolecular 
reaction and in a good qualitative agreement for bimolecular mechanism (a 
comparison of these results is provided in Appendices). Particularly, the MF 
model was able to reproduce MC results distinguishing between surfaces with 
different arrangements of active sites. Therefore, we employed the mean-field 
approximation with inhomogeneity parameters in the further calculations. 

Fig. [2] presents turn-over number as a function of surface diffusion prob- 
ability for surfaces with different arrangements of active sites in the case of 
unimolecular reaction for half of centers being active. The arrangement of 
active sites strongly influences catalytic activity of the system resulting in 
different TON values for different types of the arrangement. As the order- 
ing of active sites increases (in a series: '1+1', '2+2', '3+3', '4+4', '1/2') 
the values of TON decrease. The random arrangement leads to TON values 
located between these for '1+1' and '2+2'. For low values of Pdiff the differ- 
ences in TON are relatively small because surface diffusion is slow and the 
adsorption/desorption balance is the main factor influencing the rate of the 
overall process. As Pdiff increases the differences become more pronounced 
and '1+1' and 'rand' systems are clearly the most efficient ones. Particu- 
larly important is the case of '1/2' arrangement where TON is constant over 
the whole Pdiff range and TON values are the lowest. In each case half of 
surface centers are inactive in the reaction, however, they are active in the 
adsorption, desorption and diffusion processes. Inactive sites in the presence 
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of surface diffusion constitute an additional channel transporting reagent 
molecules from the gas phase onto active sites and thus cooperative effects 
are observed. Similar cooperative effects between active and inactive surface 
centers are observed, for example, for supported metal catalysts [3l |18] . As 
the ordering of active sites arrangement increases, the amount of active sites 
neighboring inactive ones decreases and this additional diffusive transport 
becomes less efficient. For '1/2' system this transport is negligible, since the 
length of the border between active/inactive asymptotically goes to zero. 

In Fig. [3] the values of turn-over number vs. diffusion probability for 
bimolecular mechanism are presented. Here TON is lower than for the uni- 
molecular case, however, qualitatively this plot resembles unimolecular case 
with the exception of both very low Pdiff values. For bimolecular reaction, 
the presence of pairs (at least) of neighboring active sites is a necessary con- 
dition for surface reactivity. For Pdiff (see inset) TON values for all 
systems, except 'rand', are located in the same point. For all these systems 
there exist a relatively high number of active-active neighbors whereas in the 
'rand' case this number is lower and hence the latter surface is less reactive. 
As surface diffusion becomes more significant 'rand' system is more reactive 
and above Pdiff = 0.012 it is the second reactive one. For Pdiff < 0.007 
('rand' more efficient than '1/2') TON values are determined mainly by the 
number of active-active neighbors whereas above this threshold the cooper- 
ative effects between active and inactive sites become more important and 
under these conditions active-inactive neighborhoods are also promoting the 
reaction. This cooperative effects are observed for all systems except, like 
in the unimolecular case, '1/2' one where the influence of surface diffusion is 
neglegible. 

In Fig. Ha difference between TON of 'rand' and '1/2' system (ATON) 
vs. both 9act and Pdifj is presented. This quantity depicts the difference in 
catalytic activity of the system with active sites spreaded randomly among 
inactive ones and the surface with active sites accumulated in one surface 
region. For low Pdiff there is no influence of the active sites arrangement on 
the rate of catalytic process (ATON 0). As Pdifj increases, the difference 
between 'rand' and '1/2' system becomes more pronounced and 'rand' system 
(the one with the possibility of cooperation between inactive and active sites) 
becomes more effective. There is an optimal value dact ~ 0.32 for which 
ATON achieves a maximum. At low 6act, in both 'rand' and '1/2' system, 
TON is near zero, therefore, the values of ATON are also low. While 6act 
becomes higher, the cooperative effects in 'rand' system no longer influence 
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the catalytic process (the number of inactive sites is very low) and hence 
there is no distinction between both systems. We observed similar behavior 
in the case of bimolecular reaction mechanism (results not presented here), 
in that case the maximum of ATON was located at Oact ~ 0.45. These 
results support the above conclusions about the role of surface diffusion in 
cooperative effects between active and inactive sites. 

4 Conclusions 

In this work we investigated theoretically the role of surface diffusion on the 
chemical kinetics for processes catalyzed by inhomogeneous surfaces. We 
employed extended mean-field model which takes into account the arrange- 
ment of active sites on the surface and we found analytical solutions for two 
model processes on different surfaces. We demonstrated that surface diffu- 
sion strongly influences the kinetics and that diffusion is responsible for the 
differences in the reactivity of systems with different arrangements of active 
sites. The presence of cooperative effects between inactive and active cen- 
ters was demonstrated. Inactive sites, due to the possibility of adsorption 
and surface diffusion of particles on them, constitute an additional channel 
transporting reagent particles onto active sites. The important message of 
this study is that these phenomena can be described with mean-field model- 
ing. It should be stressed that the results presented here should be treated 
quahtatively. This study shows the general features of catalytic systems with 
inhomogeneous surfaces and the possible phenomena which may be observed 
and we focused solely on the role of the arrangement of surface sites in the 
presence of finite surface diffusion. It should be noted that the mean field 
description presented here is able to reproduce qualitatively, for the consid- 
ered uni- and bimolecular processes, the whole variety of phenomena (such 
as spillover /reverse spillover, capture zone effects, support interactions and 
communication effects) that are observed experimentally for more complex 
systems (the reader can find a good summary of such experimental results 
in the review of Libuda and Freund [23] ) ■ 
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Appendices 



A Solutions of probabilities evolution equa- 
tions 

Introducing: [A*] = X, [Ao] = y, Pads,act = Pdes,act = h, Pads.inact = h, 

Pdes,inact = ^4; PdiffJA = ^5) Pdiff,AI = ^6) Pr = ^7) ^act = C-, a^mact = b 

the following set of nonlinear equations can be written for unimolecular mech- 
anism in the place of Eqs. (Cj) i ([H]): 

^ = ki (a — x) + k^wy (a — x) — 

—k2X — kQWX {h — y) — kjx (12) 

dy 

— = ks{b-y) + kewx {b - y) - 
at 

—k^y — k^wy (a — x) (13) 
For bimolecular mechanism we get: 



— = ki {a — x) + k^wy (a — x) 



dt 



—k2X — kgwx {h — y) — 2kjWrX (14) 



dy 

— = k^{h-y) + kQWX {b - y) 



dt 



—k^y — k^wy^a — x) (15) 



Here we are interested in the stationary solutions of one-site probabilities 
evolution equations (we are going to study steady-state kinetics) hence we 
assume: 

dx dy , ^ 

— = and ^ = 16) 
dt dt ^ ' 

An analytical solutions for different values of parameter w can be found 

assuming that the probability of diffusion between active and inactive region 

is equal, i.e., k^ = k^. Below are the solutions for uni- and bimolecular 

mechanism. 
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A.l Unimolecular mechanism 



X = Z y\ahk^ + a kij k^w + aki (A;4 + A^a) j (17) 
y = Z'^ {^{b"^ k3 + abki^ k^w + bkskj + bk2k3 + bkik3^ (18) 

where: 

Z = {ak^ki + {b {ki + k^) + ak2 + a ki) k^) w + 

{ki + ks) kj + k2 {ki + kz) + ki {k^ + k^) (19) 

A. 2 Bimolecular mechanism 

In the case of bimolecular reaction mechanism two solutions of the set of Eqs. 
f|T4|) and f|T5l) were obtained with opposite sign of x. The one with x > was 
taken into account (we do not present the solution for y here since it was not 
used during TON calculations): 

X = (((8 * * b * ks + 8 * * ki) * kl * k^ * + 

({8 * a * b * k^ + 16 * * ki) * k^ + 8 * a * b * kl + 16 * * ki * k-^) * k^ * ky * w + 

(8 * a * ki * kl + 16 * a * ki * k^ * k^ + 8 * a * ki * k^) * kj) * Wr + 

{b^ * k4:^ + {2*b^ * k3 + 2*a*b* k2 + 2*a*b* ki) * k4^ + b^ * kS^ + 

(2 * a * b * k2 + 2 * a * b * ki) * k^ + * kl + 2 * * ki * k2 + * kf) * kj * w"^ + 

({2 * b * k2 + 2 * b * ki) * kl + ({A* b * k2 + 4:* b * ki) * ks + 2 * a * kl + A* a * ki * k2 + 

2 * a * kl) * k^ + (2 * b * k2 + 2 * b * ki) * kl + 

(2*a*kl + 4:*a*ki*k2 + 2*a* k^) * k^) * k^ * w + 

{kl + 2 * ki * k2 + kl) * kl + {2 * kl + A* ki * k2 + 2 * kl) * ks* k4 + 

{kl + 2 * ki * k2 + kl) * klY^'^ + (—6 * k^ — b * k^ — a * k2 — a * ki) * kr, * w + 

(-/cs - ki) *k4 + (-/C2 - ki) * ks)/ 

({A * a * k^ * kj * w + {A * ki^ + A * k^) * kj) * Wr) 



B Comparison with Kinetic Monte Carlo re- 
sults 

In order to test the usefulness of the presented mean field model it was 
employed to describe the kinetics of both a uni- and bimolecular catalytic 
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process on the model surfaces with various arrangements of active centers. 
We started with the lattice models of surfaces with randomly arranged active 
sites, then the surfaces with the stripes of active sites of increasing width 
were considered. The employment of stripes let us control the ordering of 
the active centers arrangement. In each case the ratio of the number of active 
to inactive sites was kept constant {OacUve = 0.5). The same surface models 
were used in our previous Kinetic Monte Carlo study [22] • 

To make possible the direct comparison of results obtained using the 
present model with the results of Monte Carlo simulations the parameters 
of the former (the constants in Eqs. (fTTl) - (fT9l) and in the solution of Eqs. 
f ll4p and (1151) ) were chosen in a correspondence to the parameters used in 
the simulations 0. 

The crucial step in the implementation of the present model is the deriva- 
tion of inhomogeneity parameters w and Wr- The former parameter is an 
average probability that, if one considers a pair of neighboring centers, an 
active-inactive pair is found. In order to derive w the neighborhood has to 
be defined. Here we consider 8 nearest lattice sites to a given one as its 
neighborhood (this is a minimum neighborhood that distinguish among con- 
sidered arrangements of active sites because an alternative, 4 nearest lattice 
sites neighborhood, gives w that does not differ for 'rand' and '1-1-1' sys- 
tems). For the random distribution of active sites (and for OacUve = 0.5, like 
in the present case) w is equal to 0.5 because the probabihty, that an active 
center has an active neighbor, is 0.5 since they are spread randomly. For 
'1-1-1' system, w = 6/8 because among 8 neighbors of each active site 6 of 
them are inactive. For the systems: '2-1-2', '3-1-3' and '4-1-4' the values of 
w are, respectively, equal to: 3/8, 2.0/8, 1.5/8. In the two latter cases the 
neighborhoods are different for centers located inside and on the edges of 
active stripes therefore the average weighted with respect to the number of 
differently located sites was taken (it is depicted by the decimal representa- 

^In the Monte Carlo simulation algorithm different number of processes was possible for 
different lattice sites and hence the probabilities of processes were varying between active 
and inactive centers. In order to include proper values of parameters the probabilities were 
rescaled in the following way (primed values were actually put into equations whereas in the 
text the values referencing to the unprimed parameters are given): k'^ = Xgas^ ^2 = \pdes-i 
k's = Xgas, ^4 = ^pdes, K = \pdiff, K = l^Pdiff- In Order to use fcg = fcg approximation 
the values of k'^ and fcg were averaged and the arithmetic mean value fcg = O.lOApdif f was 
taken into account in calculations. We also, in order to compare TON values, renormalized 
active sites coverage [A*]' — 2 [A*] in order to get the coverage of active sites and not the 
coverage of the whole surface. 
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tion of the numerators 2.0 and 1.5). For '1/2' system in the hmit of infinite 
lattice w ^ because the number of active-inactive neighbors is neghgible 
in comparison with the total number of pairs. According to the values of w 
the considered systems can be ordered with the increasing 'inhomogeneity' 
as the following: '1/2', '4+4', '3+3', '2+2', 'random', '1+1'. It is worth to 
note that '1+1' system is , according to our definition of 'inhomogeneity', 
less ordered than the system with the random distribution of active centers. 

The derivation of the parameter Wr is less straightforward. Considering 
the bimolocular reaction, we assume that 4 nearest lattice sites form a site's 
neighborhood (the same assumption was taken in the Monte Carlo simulation 
algorithm). The rate of bimolecular surface reaction, from the point of view 
of a single active site, is proportional to the number of active centers in 
the neighborhood. Because of this it seems, that the more active neighbors 
are present, the higher the rate should be. However, active centers in the 
neighborhood are also competitors of the considered active site in the usage 
of the reactant and due to this effect the rate is, at the same time, inversely 
proportional to the number of active centers in the neighborhood. Therefore 
there exists a balance between the cooperative and competitive aspects of 
active-active neighborhood. On the other hand, at least one active neighbor 
is necessary for a given active center to be able to generate product, and 
hence 1 active site in the neighborhood is a threshold for reaction, and the 
further increasing of the number of active neighbors above 1 does not lead to 
the increase of the reaction rate. Taking into account these considerations, 
Wr can be approximated as equal to the average probability of finding at 
least one active center in the nearest neighborhood of an active site (note 
that according to this definition, Wr does not increase when the number of 
active neighbors increases above 1). The parameter Wr was equal to 1 for each 
surface model, except the one with randomly placed active centers, where (in 
the case of 4 nearest sites neighborhood): 



and in this system it was equal to 0.9375 for 6 act = 0.5. In the above equation 
the (1 — Oact)* = Ofnactive term gives the probability that each of 4 neighbors 
is inactive one, thus the 1 — (1 — OactY term corresponds to the probability 
that at least one of neighbors is active. 

Fig. [5] shows TON vs. pr parameter for the system with the unimolecular 
surface reaction mechanism. TON is calculated employing presented ana- 
lytical model with the set of model parameters corresponding to the Monte 




(20) 
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Carlo study {9act = 0.5, Xgas = 0.05, pdes = 0.01, Pdifj = 0.5). Curves for 
various arrangements of active centers, i.e., various values of w parameter, 
are presented. The results are in good agreement with the corresponding de- 
pendencies showed in [22]. The highest values of TON are reached for '1+1' 
system then the catalytic efficiency decreases along with the increasing or- 
dering in the arrangement of active sites. TON for '1/2' case is, consequently, 
the lowest one. There is a minor difference between Monte Carlo and analyt- 
ical model: in the latter one TON of 'random' system is nearer TON values 
for '2-1-2' system, whereas it was more similar to the '1-1-1' case in the Monte 
Carlo study. Also, the values of TON obtained with the presented model are 
slightly lower than those predicted in simulations. 

In Fig. [H]the dependencies of TON on pr for systems with the bimolecular 
surface reaction are presented (for the same set of parameters as in the uni- 
molecular case). These results are in qualitative agreement with our previous 
Monte Carlo study, i.e., like in the unimolecular case, TON decreases along 
with the increasing of ordering of active centers on the surface. The system 
with randomly arranged active centers is an exception, because its efficiency 
is comparable with '2+2' system, whereas TON was near '4+4' system in the 
results obtained in simulations. From the quantitative point of view, the val- 
ues of TON for the bimolecular reaction mechanism are lower in comparison 
with those obtained in simulations. The reasons of this differences between 
the mean field model and Monte Carlo simulations lay, most probably, in the 
derivation of Wr parameter wherein the size of considered neighborhood can 
be chosen in different ways. Since in this work we are mainly interested in 
the qualitative description of considered catalytic systems, we may conclude 
that in both unimolecular and bimolecular case the mean field model is able 
to reproduce the most important features of deliberated systems. 
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Figure 1: Schematic representations of considered surfaces, white and gray 
lattice sites represent, respectively, inactive and active surface centers. 
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Figure 3: Turn-over number vs. probability of surface diffusion for bimolec- 
ular reaction idact = 0.5, Xgas = 0.05, pdes = 0.01, pr = 0.9). 





Figure 4: ATON vs. fraction of active sites and surface diffusion rate for tlie 
unimolecuar surface reaction mechanism {xgas = 0.05, pdes = 0.01, Pr = 0.9). 
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Figure 5: TON vs. reaction probability for the unimolecuar surface reaction 
mechanism. 
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Figure 6: TON vs. reaction probability for the bimolecular surface reaction 
mechanism. 
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